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The soft behavior of the bremsstrahlung from a source is discussed in terms of classi- 
cal transport models and within a non-equilibrium quantum field theory (Schwinger 
- Kadanoff - Baym - Keldysh) formulation. 



Introduction 

(N . 

We study the importance of coherence time effects (Landau - Pomeranchuk - Migdal - ef- 
(*C) 1 feet (LPM)|[|l|) on the production and absorption of field quanta from the motion of source 

particles in non-equilibrium dense matter. In order to calculate such bremsstrahlung ef- 
fects appropriately one needs to go beyond the commonly used quasi-particle picture, 
and include effects of the finite widths of the particles fj, [|, £§]. Our considerations are 
of particular interest for the application to photon, or di-lepton production in high en- 
ergy nuclear collisions, for gluon or parton radiation and absorption in QCD transport 
and its practical implementation in parton kinetic models (such problems are discussed, 
e.g. in H, |5[), to neutrino and axion radiation from supernovas and neutron-star matter 
(see ||), for the soft phenomena in quantum cosmological gravity (see 0) and also for 
many condensed matter phenomena, as particle transport in metals and semiconductors, 
radiation in plasma etc (see ||). The study also gives some hints how to generalize the 
standard transport picture (see ||) such that it can include genuine off-shell effects due 
to damping of the single particle propagation from collisions or decays. 

To be specific we take the example of electrodynamics, considering photon production 
and assume that the source system couples only perturbatively (to lowest order) to the 
radiation field, while the source itself can interact in any non-perturbative way. Thus, 
we consider a "white" body as a source for the radiated field! Our considerations are 
formulated in real time non-equilibrium Green's function technique, where the production 
rate is given by the — h component of the proper self-energy diagram of the produced 
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photon 

- iIT + = ^(Jjll)-- = Arc I d A ie^ (f\x - i/2)f(x + £/2)) , (1) 

which is determined by the current auto correlation function of the source. The dashed 
lines relate to the photon, while the —ill-loop symbolically denotes the exact inclusion of 
all strong interaction among the source particles. The bracket (...) denotes a quantum 
ensemble average over the source with quantum states and operators in the interaction 



picture. Above and throughout later we use the Keldysh {— , +} notation, see |10fl , where 
the n~ + and IT + ~ self energies are responsible for gain and loss. Such a formalism has 
been applied in calculation of neutrino emissivity of neutron stars in ref. [nj, employing 



the quasi-particle approximation for the equilibrium nucleon Green functions. However 
the general formalism allows to go beyond this limit and to account for the finite damping 
width of the source particles due to their finite mean free path. 

Classical radiation 

For some classical source systems perturbatively coupled to a boson field the boson self 
energy can be obtained in closed form. For these cases one has to evaluate the current- 
current correlator on the classical level. The properties of classical radiation will be 
illustrated for two examples, where a charged particle (the source) stochastically moves 
in dense matter. Thereby the motion of the charge is described either (a) by mesoscopic 
transport (diffusion process) or (b) by a microscopic Langevin process. 

Diffusion process 

The motion of the source particle is assumed to be described by a time dependent 
phase-space distribution f(x,v,t) in space and velocity with convective current density 
j{x,t) = e J d 3 v v f(x,v,t). For standard dissipative media in equilibrium the velocity 
autocorrelation function (integrated over space) decays exponentially in 

fc^N = I /,7 2 \ wk e -r x \r\ 



(r)v k (0)) = -(v 2 )5 
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where i and j denote the spatial components. T x is the relaxation rate which we approx- 
imate as constant, (...) now denotes the average over the classical distribution functions. 

Both f(x,v,t) and the autocorretion function can be obtained in closed form, if the 
time evolution of /, and the propagation of fluctuations 5f are governed by a standard 
diffusion process. Solving then a Fokker-Planck equation for the space-time dependence 
of velocity fluctuations we obtain (in mixed r, q representation) 



-m;+(r,q) = Arc | d^e"^ (/(f , r)/(0, 0)^ 

= Arce 2 Po {((v*v k )) eq e" r ^ - D^q k ( e - r -M - l) 2 } (3) 
x«p|-^(r,|r|+e- r -W-l)}. 



4 We should note that this ansatz ignores finite size corrections in terms of some anti-correlations on 
the scale of a mean recurrence time Hi . 
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Here the ensemble average (. . .) eq over the equilibrium distribution f eq keeps only even 
moments of v with (v 2 ) = 3DT X such that only even powers of iq survive. This renders 
— iil^ + (r, q) real and symmetric in r. For transverse photons fjq = some terms drop 
and we have 

-in- + (^, q) = Ane 2 p (vV) eg exp [D <f/T x ] 

x f^j ( -D? \ k 2(k + l)T x + 2D<f 

^k\ \ T x J ^ k + l)T x + Dq 2 f + uj 2 
— Ane 2 Po (v^) eq 2T x /(u 2 + Tl) for (^<I^. (5) 

Compared to the infra-red divergent quasi-free result oc 1/uo 2 this form of the correlation 
renders the photon self energy regular in the soft limit at four momentum q = 0. It is 
determined by mesoscopic transport properties, namely by the relaxation rate T x and the 
diffusion coefficient D = (v 2 ) /3T X ; po is the spatial density of the charged particles. For 
large q, i.e. Dq 2 ^> T x , only the short time behavior of the autocorrelation matters and 
from (^) and (Q) one finds in (r, q) and (u, q) representations 

ton [-iU-+(T,q)] = 47re 2 P o {v l v k ) exp [-T x \r\ - D fT x r 2 /2] , (6) 



3/2 < ^ 



( 2tt \ 1 

D lim_[-in-( W , ? -)] = tae-nW^ ) ex Pl ^ 



(7) 



Microscopic Langevin process 

In a microscopic Langevin process hard scatterings occur at random with a constant 
mean collision rate T. These scatterings consecutively change the velocity of a point 
charge from v m to v m+ i,. . . (subscripts m, and n below refer to the collision sequence). In 
between scatterings the point charge moves freely. For such a multiple collision process 
one can at least determine the q = part of the self energy. 

The modulus of the autocorrelation function takes a Poissonian form for such a colli- 
sion sequence 

°° \Yr\ n 

-iU-+(r,q = 0) =A7re 2 p (v\r)v k (0)) = A^e 2 p^^^ {<<+n) m , (8) 

n=0 

where (. . .) m denotes the average over the discrete collision sequence {m}. The time 
Wigner transform of (|8|) determines the spectrum at vanishing q for all u 

n . . 2 X 2 ^/ i i \ 2r n Re{(r + icj) n+1 } 
- d£>, q = 0) = Arre 2 p £ <« + n> m V (9) 

n=0 ^ ' 

This is a genuine multiple collision description for the photon production rate in com- 
pletely regular terms due to the (u 2 + T 2 ) n form of all denominators. Each term is regular, 
since right from the beginning one accounts for the damping of the source particle be- 
cause of the finite mean time 1/T between collisions. The result (|9|) still accounts for the 
coherence of the photon field, now expressed through the correlations (v m v m+n ) m arizing 
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from the sequence of collisions. The terms in (pj) define partial rates, which are associated 
with specific self energy diagrams. 

For the Langevin process the ^-dependence of the self energy cannot be obtained in 
closed form in general. Still the n = term (c.f. with n = term from eq.@) can be 
given 

-mrM'^ ^^p ^ do) 

It shows the typical Cherenkov enhancement at oj — qv. Although for q — > the analytical 
form of ( |10D resembles the diffusion result @, it is not the same unless (v m v m+n ) m = for 
n 7^ 0, an approximation recently used in ref. [JTJ]. In the general case velocity correlations 
between successive scatterings exist, and there will be a sizeable difference between the mi- 
croscopic mean collision rate T and the mesoscopic relaxation rate T x . For systems, where 
the velocity is degraded by a constant fraction a, such that (v m ■ v m+n ) m = a n (v m ■ v m ) m , 
one can sum up the whole series in (BJ) and recover the diffusion result (PS) at q — with 
T x = (1 — a)T. This clarifies that the diffusion result (|]) represents a resummation of the 
random multiple collision result @. 

Monte Carlo evaluation of amplitudes 

Some of the cascade schemes try to cure the infra-red problem by considering the phase 
of the photon field along the classical orbits. Thus they evaluate f dtv(t) exp[io;t — iqx(t)} 
along the random straight sections of the classical paths in the cascade model. We like to 
mention that this is a highly unreliable procedure due to the strong cancelations of terms 
that are randomly generated. Our experience shows that a reliable result (better than 
< 5%) in the transition region u w V requires an ensemble of 10 3 cascade runs where each 
path has about 10 3 collisions. Compared to that the analytical result (||) has significant 
computational advantages, provided the random process is of this form. 



Quantum many-body description 

In the full quantum field theory formulation the production rate is given by all photon 



self energy diagrams of skeleton type according to the closed time-path rules |H| [L3[ with 
full Green's function which result from the summation of Dyson's equation. The later in 
short matrix notation becomes 

G = G O + G O 0S0G (11) 



+ -iS 



Here G , G and S denote the two by two matrices of the unperturbed Green's function 
(thin lines), the full Green's function (full lines), and the proper self energy of the source 
particles. The denotes the space-time folding. The four components of — iS are defined 
as the sum of all standard proper self energy diagrams like in normal perturbation theory, 
now however with definite + or — assignments at the external vertices, and summed over 
the — and + signs for all internal vertices. These signs specify the Green's functions 
linking the vertices and the choice between the — vertex — iV and its adjoint value +iV^ 
at + vertices (c.f. ref. El). 
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The full Green's functions account for the finite damping width of the particles, and 
therefore destroy of the strict energy-momentum relation in dense matter. This width 
expressed through the imaginary part — r of the retarded self energy arises from collisions 
or decays of the particles. In this picture the off-diagonal Green's functions G~ + and 
G + ~ are four momentum and space Wigner densities for the occupied and available 'single 
particle states', which now have a finite width. The corresponding off-diagonal parts in 
the self energy determine the corresponding gain and loss terms in a transport description. 
Therefore each diagram of ([J) defines a specific partial rate. 

We suggest a decomposition of the diagrams that allow for a simple interpretation and 
classification in terms of physical in-medium scattering processes, and propose particular 
resummations of physically meaningful diagrams, which consider the finite damping width 
of all source particles in matter. All the graphs consisting G~ + G + ~-products are explicitly 
presented. In this picture the set of diagrams reduces to 



+••• (12) 



Here full dots and boxes denote effective in-medium vertices and 4-point interactions, e.g. 

J | = j + j | | + jnf ' ••• (is) 

These in-medium vertices and four point interactions are defined for a specific choice of 
sign (say — ) through skeleton resummation schemes, where only bare — vertices linked by 
G Green's functions appear. Compared to conventional diagrams, vertex corrections 
of different signs appear on both sides of a loop once they are separated by {H — } lines. 
In ( |i~3"D the wavy lines are either interactions or full boson-propagators in a theory of 
fermions interacting with bosons, like in QCD. In some simplified representations (being 
often used, c.f. ||11|| ) the 4-point functions behave like intermediate bosons (e.g. phonons). 

We note that each diagram in (|T^) represents already a whole class of perturbative 
diagrams of any order in the interaction and in the number of loops. The most essential 
term is the one-loop diagram in (0), which is positive definite, and corresponds to the 
first term of the classical Langevin result for Il c i in (|9|). The other diagrams represent 
interference terms due to rescattering. All diagrams calculated with full Green's functions 
are void of infra-red divergences. Thus these diagrams represent the quantum general- 
ization of the infra-red regular Langevin result ([|). Under appropriate conditions, the 
correct quasi-particle (QP) and quasi-classical (QC) limits are recovered from this subset 
of graphs. 

Decomposition of closed diagrams into Feynman amplitudes in the quasi-particle approx- 
imation (QPA) 

The QPA is a quite commonly used concept originally derived for Fermi liquids at low 



temperature (Landau - Migdal, see JL4]), where it constitutes a consistent scheme. With 
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the application of transport models to high energy situations this concept has been taken 
over to a regime where its validity cannot be justified under all circumstances. For the 
validity of the QPA one normally assumes that T e, where e is an average particle kinetic 
energy (~ T for equilibrium matter). This has considerable computational advantages as 
Wigner densities (" — h" and "H — " lines) become energy (^-functions, and the particle 
occupations can be considered to depend on momentum only rather than on the energy 
variable. |. Formally the energy integrals can be eliminated, in diagrammatic terms just 
cutting the corresponding " — h" and "H — " lines [|ll| . 



Thus the QP picture allows a transparent interpretation of closed diagrams. With 
consecutive numbers 1 to 6 for the diagrams drawn in ([12]), diagrams 1, 2, 4 and 5 relate 
to the radiation from a single in-medium scattering between two fermionic quasi-particles 
jnJ . Thereby diagram 2 is more important than diagram 4 for neutral interactions, while 
this behavior reverses for charge exchange interactions (the latter is also important for 
gluon radiation in QCD transport due to color exchange interactions). Diagrams like 3 
describe the interference terms due to further rescatterings of the source fermion with 
others. For diagram 6 the photon is produced from intermediate states (its contribution 
is suppressed in the soft limit). Some of the diagrams, which are not presented explicitly 
in eq . (|T2"1) give more than two pieces, if cut, so they do not reduce to the Feynman 
amplitudes. Since one works with zero-width fermion Green's functions in QPA, the 
finite width contributions have to appear in higher order diagrams through corresponding 
ImE-insertions! Therefore the whole set of diagrams defining the full — in h in QPA is 
by far larger, than ours (|T2"D. 

Moreover, our considerations show that the validity condition T e is not sufficient 
for the QPA. Rather, since finally energy differences of order uj appear, one has to demand 
that also uj 3> T in the QPA. In particular, the remaining series of QP-diagrams is no 
longer convergent unless uj > T, since arbitrary powers in V /uj appear, and there is no 
hope to ever recover a reliable result by a finite number of QP-diagrams for the produc- 
tion of soft quanta! With full Green's functions, however, one obtains a description that 
uniformly covers both, the soft (uj <C T) and the hard (r <C uj) regime. 



Quasi-classical limit (QCL) 

The QCL requires that i) the particle occupations are small ({np) <C 1) implying a 
Boltzmann gas and that ii) all inverse length or time-scales times h are small compared 
to the typical momentum and energy scales of the source systems. In particular we shall 
assume fauj <C e, and a fermion collision rate T = h/r co n <C e. The last inequality allows to 
use Kadanoff-Baym ansatz, where the fermion occupations depend only on momentum 
rip- = n tfl -^ F but no longer on the energy e of the particles. Assuming again V ~ const. 
on relevant time scales one immediately finds 

G~ + (r,p) ~ irvexp[-|rr|/2 + ie^r], (14) 

G+-(r,p) ~ -i(l - n ? ) exp[-|rr|/2 - ie^r]. (15) 
With these Green functions we calculate the diagrams (|T2|). For the one-loop diagram we 

5 The later approximation is also often used beyond the scope of the QPA and is then known as 
Kadanoff-Baym ansatz (Tj|, c.f.]8| 
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obtain the expression identical to the n = term of the classical Langevin result (|T0|): 



° ^~ 4?re 7 ( 27 r)3 (^-^f + P 47reP0 \(^-gl;)) 2 + rV ' W 

For neutral interactions (corresponding to the classical examples) we show that precisely 
the diagrams in the first line of (|T^) denoted by — ill~ + with n { — +} loop insertions 
correspond to the n-th term of the Langevin result (M)f\ Thus 



irvi n 

-in- + (r,g = 0) » 47re 2 p (vlv k m+n )J—]-e-^\, (17) 
-iLW,g = 0) « 4vre 2 po (<^ + n) m 2r"Re _ ^ . (18) 



For neutral interactions all other diagrams in ( |T2|) drop in the QC limit. Diagram 
4 is unimportant for neutral interactions. Other diagrams in the original series (|TJ) 
acquire extra powers either in the mean occupation (np) C 1 or in T/e < 1. Such 
suppression factors result from extra G h lines compared to the classical diagrams of 
the same topology or are due to a violation of time ordering (diagram 5), since classical 
interactions are of time-scale 1/e, which is short compared to the damping time 1/T X . 

We note that truncating the set of diagrams (|12|) with respect to a certain expansion 
parameter special care should be taken in order to satisfy charge-current conservation 
law. For verification one may use the Ward-Takahashi identities within the desired order. 

For equilibrium T ^ case even in general quantum consideration the one-loop dia- 
gram can be expressed by the QP prescription [JIT], [151 however multiplied by the factor 
C = uj 2 /{uj 2 + T 2 ), which displays the suppression at low u. There is hope that even in 
the quantum case some higher order diagrams can also be resummed and that an over- 
all suppression factor of the form C = uj 2 /(uj 2 + T 2 ) emerges for the true in-matter rate 
relative to the quasi-free or QP one in the limit q = 0, c.f. the diffusion result (|[-^). 

Altogether our results provide an extension of the incoherent quasi-free and QP ap- 
proximation from the ultra violet limit towards the soft limit with appropriately modified 
production cross sections. It does not only regularize the infra-red divergence of the free 
rate, but it also produces the right q dependence in the soft limit. We also note that 
our diagrammatic description may suggest a formulation of a transport theory which in- 
cludes the propagation of particles with finite width and therefore may permit a consistent 
treatment of resonances in non-equilibrium dense matter. 



6 The proof is easily given in the r — p representation, assigning times and r to the external — and 
+ vertices, while the internal — and + vertices are taken at t± to t~ and t\ to i+, respectively. In the 
classical limit G is retarded, while G ++ is advanced, such that both time sequences have the same 
ordering: 0<i^<...<<^<r and < <...<<+< r (the order reverses, if the line sense of 
the outer fermion lines are reversed). Thus the r-dependence of the modulus of these diagrams again 
results in e~' rT L Assuming classical momentum transfers \p n — p n +i\ which are large compared to HT one 
concludes that the {H — } loop insertions merge S(t~ — f+) in their time structure. With t n = t~ — t+ a 
given diagram then no longer depends on the intermediate times t n apart from the ordering condition, 
and therefore results in a factor |rY|™/?i!. With q = also the corresponding momenta are pair-wise 
identical, and the remaining momentum integrations just serve to define the correlation between v m and 
v m +n after n collisions of the charged particle. 
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